Fluid-solid phase transitions in 3D complex plasmas under microgravity conditions 
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Phase behavior of large three-dimensional complex plasma systems under microgravity conditions 
onboard the International Space Station is investigated. The neutral gas pressure is used as a control 
parameter to trigger phase changes. Detailed analysis of structural properties and evaluation of three 
different melting/freezing indicators reveal that complex plasmas can exhibit melting by increasing 
the gas pressure. Theoretical estimates of complex plasma parameters allow us to identify main 
factors responsible for the observed behavior. The location of phase states of the investigated 
systems on a relevant equilibrium phase diagram is estimated. Important differences between the 
melting process of 3D complex plasmas under microgravity conditions and that of flat 2D complex 
plasma crystals in ground based experiments are discussed. 
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INTRODUCTION 



Complex (dusty) plasmas - systems consisting of 
highly charged micron-size particles in a neutralizing 
plasma background - exhibit an extremely rich variety of 
interesting phenomena [l|-|8| . Amongst these transitions 
between fluid and solidphases (freezing and melting) are 
of particular interest [l|, y, |j, 



gas outlet 



14|. This is largely a 



consequence of the fact that high temporal and spatial 
resolution allows us to investigate these phase changes 
along with various related phenomena at the individual 
particle level [1, 4, 7, 8, 13 22J. Fully resolved atomistic 
(undamped) particle dynamics provides new insight into 
natural atomic and molecular systems, whose dynamics 
cannot be resolved in such detail. 

In this paper we report experimental investigations 
of the fluid-solid phase transitions in large 3D com- 
plex plasmas performed under microgravity conditions 
onboard the International Space Station (ISS). These 
phase changes are driven by manipulating the neutral 
gas pressure. Detailed analysis of complex plasma struc- 
tural properties allows us to quantify the extent of or- 
dering and accurately determine the phase state of the 
system. Evaluation of various freezing and melting in- 
dicators gives further confidence regarding the phase 
state. It is observed that the system of charged parti- 
cles can exhibit melting upon increasing the gas pres- 
sure, in contrast to the situation in ground-based exper- 
iments where plasma crystals normally melt upon reduc- 
ing the pressure. This illustrates important differences 
between generic (e.g. similar to conventional substances) 
and plasma-specific mechanisms of phase transitions in 
complex plasmas. 

First results from the studies described here have been 
reported in Refs. 23l l2J|. The purpose of this paper is 
to provide more detailed and comprehensive information 
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FIG. 1: (Color online) Sketch of the PK 3 Plus plasma cham- 
ber [Hi. 



on the experimental procedure, analysis of the obtained 
results, and their theoretical interpretation. 



EXPERIMENTAL SETUP 

Experiments are performed in the PK-3 Plus labora- 
tory onboard the ISS [25|]. The heart of this laboratory is 
a parallel-plate radio-frequency (rf) discharge operating 
at a frequency of 13.56 MHz, sketched in Fig. [TJ The 
electrodes are circular plates with a diameter 6 cm made 
of aluminium. The distance between the electrodes is 
3 cm. The electrodes are surrounded by a 1.5 cm wide 
ground shield including three microparticle dispensers on 
each side. The dispensers are magnetically driven pistons 
filled with monodisperse particles of various size and ma- 
terial [25|. 

Discharge can operate in argon, neon or their mixture 
in a wide range of pressures, rf-amplitudes and rf-powers. 
The working pressures are in the range between 5 and 
255 Pa with and without gas flow. The latter can be 



produced by a specially designed system, allowing us to 
operate with the lowest flow rates necessary to provide 
clean experimental conditions. The flow is essentially 
a symmetrical gas curtain around the electrode system 
flowing from the lower to the upper part (see Fig. [T|). 

Complex plasmas are formed by injecting monodis- 
perse micron-size particles into the discharge. The op- 
tical particle detection system consists of a laser illumi- 
nation system and three video cameras (there is also a 
fourth camera of the same type which is used to observe 
glow characteristics of the plasma). Two diode lasers 
(A — 686 nm) are collimated by a system of several lenses 
producing a laser sheet perpendicular to the electrodes 
with different opening angles and focal points. Three pro- 
gressive CCD-cameras detect the reflected light at 90°. 
An overview camera has a field of view (FoV) of about 
60 X 43 mm^ and observes the entire field between the 
electrodes. A second camera has a FoV of about 36 x 26 
mm^ and observes the left part of the interelectrode space 
(about half of the entire system) . The third camera is the 
high resolution camera with a FoV of about 8x6 mm^ . 
It can be moved along the central axis in the vertical 
direction. The cameras and lasers are mounted on a hor- 
izontal translation stage allowing a depth scan through, 
and, therefore, 3D observation of complex plasma clouds. 

Further details on the PK-3 Plus project can be found 
in a comprehensive review |25|. 
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EXPERIMENTAL PROCEDURE 

The experiments described here are carried out in 
argon at a low rf-power (~ 0.5 W). We use two dif- 
ferent sorts of particles in the two distinct experimen- 
tal runs: Si02 spheres with a diameter 2a — 1.55 /im 
and Melamine-Formaldehyde spheres with a diameter 
2a = 2.55 /im. The experimental procedure, identical 
in these two runs, is as follows: When the particles form 
a stable cloud in the bulk plasma, the solenoid valve to 
the vacuum pump is opened, which results in a slow de- 
crease of the gas pressure p. Then, the valve is closed and 
the pressure slowly increases due to the gas streaming in. 
(Neutral flow has negligible direct effect on the particles). 
During the pressure manipulation (~ 6 minutes in total), 
the structure of the particle cloud is observed. The ob- 
servations cover the pressure range from p ~ 15 Pa, down 
to the lowest pressure of p ~ 11 Pa and then up to p ~ 21 
Pa [see Fig.[2;a)]. 

In order to get three-dimensional particle coordinates, 
30 scans are performed. Scanning is implemented by si- 
multaneously moving laser and cameras in the direction 
perpendicular to the field of view with the velocity 0.6 
mm/s. Each scan takes ~ 8 s, resulting in the scan- 
ning depth of ~ 4.8 mm; the interval between consecutive 
scans is ~ 4 s. The particle positions are then identified 
by tomographic reconstruction of the 3D-pictures taken 



FIG. 2: (Color) (a) Side view of the particle cloud (inverted 
colors) taken with the overview camera (left) and the corre- 
sponding FoV of the high resolution camera (right) [particles 
are color-coded to see solid-like (red) and liquid-like (blue) 
domains]. Rectangle marks the part of the cloud used for the 
detailed structural analysis (rectangular box 7.0 x 0.7 x 4.5 
mm^). (b) Thickness of the particle cloud in the vertical 
direction vs. the scan number. Blue triangles (red circles) 
connected by lines correspond to the system of small (large) 
particles. The corresponding values of pressure are shown by 
a brown solid curve (the dependence of pressure on the scan 
number is almost identical in the two runs). The plateaus 
on the width-curves seen when the pressure increases (right 
side of the figure) correspond to the situation when the upper 
cloud boundary leaves the FoV of the high resolution camera; 
the actual width of the cloud is somewhat larger than shown. 
(c) Mean interparticle separation A (in the part of the cloud 
chosen for the analysis) vs. the scan number in the two ex- 
perimental runs. Blue (red) color corresponds to the system 
of small (large) particles. Insets show the dependence A(p) 
demonstrating some hysteresis, which is more pronounced for 
small particles. 



with the high resolution camera observing a region 8x6 
mm" slightly above the discharge center. 



ANALYSIS 



Global reaction 



Let us first analyze the global reaction of the particle 
cloud on the pressure manipulation. An example of the 
particle cloud as seen by the overview and high resolu- 
tion cameras is shown in Fig. [DJa). Figure [IJb) shows 
the cloud thickness in the vertical direction as a func- 
tion of the scan number (time) for both the systems of 
small and large particles. It is observed that the position 
of the upper boundary is strongly correlated with pres- 
sure: It moves downwards (upwards) with the decrease 
(increase) of p. This has a clear physical explanation. 
Particles cannot penetrate in the region of strong elec- 
tric field (sheath) established near the upper electrode. 
The position of the upper cloud boundary is thus set by 
the sheath edge. The sheath thickness is roughly propor- 
tional to the electron Debye radius Adb which exhibits 
the following approximate scaling Adb oc rie oc p~^'^, 
where ne is the electron density. This implies that upon a 
decrease in the pressure, the particles are pushed farther 
to the electrode and vice versa, in full agreement with 
the observations. 

The lower cloud boundary, associated with the pres- 
ence of the particle- free region (void) in the central area 
of the discharge |26l43Cl| . shows less systematic behavior. 
Its position does not change when the pressure decreases, 
but then moves slightly upwards when it increases. How- 
ever, the displacement amplitude is relatively small. As a 
result, the thickness of the cloud exhibits pronounced de- 
crease (increase) when the pressure decreases (increases), 
as can be seen in Fig. EKb). The plateaus on the cloud 
width-curves seen when the pressure increases are arti- 
facts due to the upper cloud boundaries leaving the FoV 
of the high resolution camera (so that the actual width 
of the clouds is somewhat larger than shown). Thus, 
the particle component becomes compressed by reducing 
the pressure and expands when the pressure is increased. 
The resulting dependence of the mean interparticle dis- 
tance (in the part of the particle cloud subject to de- 
tailed analysis) on the scan number/pressure is shown in 
Fig. [2](c). The mean interparticle distance A is clearly 
correlated with pressure, although some hysteresis (more 
pronounced for the system of small particles) is evident 
from insets in Fig. ^c) . 



Structural properties 

The observed clouds of particles are not very homo- 
geneous. For example, typical interparticle separations 
in peripheral regions close to the cloud boundaries can 
exceed those in the central part of the cloud by a factor 
of about two. For this reason, a relatively small central 



part of the cloud sketched in Fig. [2ja) has been chosen 
for the detailed analysis of the structural properties. This 
part is sufficiently small (especially its vertical extent) so 
that the system inside is reasonably homogeneous. At 
the same time, it contains enough particles (~ 10'') to 
yield reasonable statistics. 

To determine the local structural properties of the 
three-dimensional particle system we use the bond or- 
der parameter method SlJ, |32|, which has been widely 
used to characterize order in simple fluids, solids and 
glasses 3l|-l34|. hard-sphere (HS) systems [35| - |4C| . col- 
loidal suspensions 4ll, 42 1, and, more recently, 3D com- 



plex plasmas 171, 120, |43|-|45| as well as complex plasma 
films 46|, |47J. In this method, the rotational invariants 
of rank I of both second qi(i) and third w/(z) order are 
calculated for each particle i in the system from the vec- 
tors (bonds) connecting its center with the centers of the 
Nnn{i) nearest neighboring particles: 
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where qim{i) = ^nn(«)~^ Ejl^i*''' ^m(i"ii), Yim are the 
spherical harmonics and r^j = r^ — r j are vectors connect- 
ing centers of particles i and ?'. In Eq.(l2l) 
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denote the Wigner 3j-symbols, and the summation in 
the latter expression is performed over all the indexes 
TTii = —l,...,l satisfying the condition mi + 7712 + "ms = 0. 
The calculated rotational invariants qi, Wi are then com- 
pared with those for ideal lattices 



4J]. Here, 



we are specifically interested in identifying face-centered 
cubic (fee), body-centered cubic (bee), hexagonal close- 
packed (hep), and icosahedral (ico) lattice types and, 
therefore, use the invariants (74, gg, wq calculated using 
the fixed numbers of N^^-^ — 12 and iV„n = 8 nearest 
neighbors, respectively. 

An example from this analysis is presented in Fig. |31 
which shows representative particle distributions on the 
plane ((74, qe). Initially, the system of small particles re- 
veals weakly ordered fluid-like structure [Fig. [31Ja)]. The 
system of large particle demonstrates more order and is 
apparently closer to the solid state [Fig. [31Jd)]. Upon a 
decrease in pressure, the particles tend to form more or- 
dered structures. At the minimum pressure both systems 
are in the solid state as can be evidenced from Figs, ^h) 
and (e) . Clear crystalline structures which are dominated 
by the hep and fee lattices are observed. Subsequent in- 
crease in the pressure suppresses the particle ordering. 
Figures [3] (c) and (f ) demonstrate the final states of the 
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FIG. 3: (Color) Variation of the structural properties with pressure as reflected by particle distributions on the plane of 
rotational invariants {q4,q(i) for the system of small (top panel) and large (bottom panel) particles. The rotational invariants 
for perfect hep and fee lattices and liquid-like domain [sketched in (a)] are also indicated. For discussion, see the text. 



systems, which are considerably less ordered than those 
in Figs. [3i;b) and (e). 



Evolution of the structural composition of complex 
plasmas in response to pressure manipulation is shown in 
Fig. 01 We see that the solid phase is mostly composed 
of hep- and fcc-like particles with only a small portion 
of bcc-like clusters. A decrease (increase) in pressure en- 
hances (suppresses) ordering of the particles. Maximum 
number of particles in the crystalline state clearly cor- 
responds to the pressure range near the minimum. The 
system of small particles exhibits melting with increasing 
pressure, which is reflected by a significant drop in the 
number of crystalline particles. No such drop is evident 
for the system of large particles, indicating that it likely 
remains in the solid phase. In this respect we also men- 
tion that a premelting stage (disappearance of fee- and 
bcc-like particles, predicted in [1J,|4J]) is observed in the 
system of small particles when the pressure increases, but 
is absent in the system of large particles. 



Freezing and melting indicators 

Several approximate approaches have been proposed 
to locate fluid-solid coexistence of various substances. 
This includes well known phenomenological criteria for 
freezing and melting, like e.g. the Lindemann melt- 
ing law 48 1, Hansen- Verlet freezing rul e 1491 . Raveche- 
Mountain-Streett criterion for freezing [50|, and a dy- 
namical criterion for freezing in colloidal suspensions 51 1 
(for a review see Ref. '52l|). These criteria are typically 
based on the properties of only one of the two coexist- 
ing phases and predict quasi-universal values of certain 
structural or dynamical quantities at the phase transi- 
tion. Quasi-universality in this context means that a 
quantity is not exactly constant, but varies in a suffi- 
ciently narrow range for a broad variety of physical sys- 
tems. It is instructive to consider application of some 
of these criteria to complex plasmas investigated in the 
present experiment. 

The Lindemann melting rule states that a crystalline 
solid melts when the root-mean-square displacement of 
particles about their equilibrium lattice positions exceeds 
a certain fraction of the characteristic nearest neighbor 




neighbor distance from its (local) mean value 
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FIG. 4: (Color online) Structural composition of complex 
plasmas composed of small (a) and large (b) particles versus 
the scan number. Relative number densities of hep-like, fcc- 
like, bcc-like and liquid-like particles (indicated in the figure) 
are shown, revealing relative abundance of different phases at 
different pressures (pressure is shown by the orange circles on 
the right-hand axis). 



distance. The critical fraction, known as the Lindemann 
parameter L, is expected to be quasiuniversal L ^ 0.15. 
In fact, however, its exact value may depend on such fac- 
tors as crystalline structure and nature (shape) of the 
interparticle interactions. For instance, for the inverse- 
power-law (IPL) family of potentials \U{r) oc r~"] the 
values of L at melting have been found to lie in the 
range between ~ 0.12 and ~ 0.15 for a wide range of n, 
3 ;$ n- < 100 53]. Somewhat higher values of L between 
~ 0.15 and ~ 0.18 have been recently reported for the 
IPL potentials with 6 < n < 10, as well as for the model 
Gaussian and exp— 6 potentials [5J]. The conventional 
Lindemann ratio is determined for the solid phase only. 
Generalizations of the Lindemann ratio which can be ap- 
plied in both the solid and fluid phases have been dis- 
cussed in Ref. [55| (see also references therein). Here we 
define the following Lindemann-like measure convenient 
for the present analysis. For each particle i in the sys- 
tem we calculate the local square deviation of the nearest 
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where dnn(i) = N^^n{i)J2f=l^'^ nj and nj = \r, - rj\. 
The global measure is then obtained by averaging over 
all the N particles in the system. 
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Since the solid phase in our case is dominated by fee and 
hep lattices we find it convenient to use a fixed number 
of nearest neighbors, iVnn — 12. The behavior of the 
Lindemann measure defined by Eq. Q upon pressure 
variations is shown in Fig. Elja). For the system of large 
particles, L is almost independent of the pressure and 
remains in the relatively narrow range L ~ 0.09 — 0.10. 
For the system of small particles, L demonstrates similar 
behavior (L ~ 0.10 — 0.12) for the first half of the obser- 
vation sequence, but then increases considerably (up to 
L ~ 0.15 at the maximum) in its second half (which cor- 
responds to an increase in the neutral gas pressure) . This 
increase in L can be interpreted as a signal of melting in 
the system of small particles. On the other hand, almost 
constant values of L for the system of large particles is an 
indication that the system remains in the solid phase. In 
addition, the fact that the values of L are systematically 
smaller for the system of large particles reveals better 
ordering in this system. All these observations are in 
full qualitative agreement with the results of structural 
analysis described in the previous section. 

The Raveche-Mountain-Streett criterion of freez- 
ing [SOf is based on the properties of the radial distri- 
bution function (RDF) g{r) in the fluid phase. It states 
that near freezing, the ratio of the values of g{r) corre- 
sponding to its first nonzero minimum and to the first 
maximum. 



R = 5(r-min)/5(r'max), 



(5) 



is constant, _R ~ 0.2 [50|. This criterion describes fairly 
well freezing of the classical Lennard- Jones fluid, but is 
also not universal. For example, in studying fluid-solid 
coexistence of the IPL systems, the ratio R at freezing 
was documented to vary between ~ 0.15 and ~ 0.25 for 
n in the range 3 < n < 100 [5^. Similarly, the values of 
R in the range between ~ 0.18 and ~ 0.24 at freezing of 
the IPL, Gaussian core and exp —6 model potentials have 
been reported 5J]. Figure. [5jb) shows the calculated val- 
ues of the freezing indicator R for different scans. Ap- 
plying the threshold condition i? ~ 0.2 would imply that 
the system of small particles melts upon an increase in 
the neutral gas pressure (second half of the observation 
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FIG. 5: (Color online) Comparison of various freezing and 
melting indicators for complex plasmas composed of small 
(blue triangles) and large (red circles) particles. The Lin- 
demann measure (a), the Raveche-Mountain-Streett ratio (b) 
and the melting indicator M (c) defined in the text are shown 
for each scan number. For discussion, see the text. 



sequence) , while the system of large particles remains in 
the solid state. This is consistent with the results from 
previous analysis. 

It has been shown recently that the cumulative distri- 
bution function 
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is very sensitive to the phase state of various systems [IJ, 
|40|. In equation ([6]), n{wQ) is the distribution, normalized 
to unity, of particles over the rotational invariant wq. 
Figure [S] shows n{wQ) and the corresponding We{wQ) for 
four different scans corresponding to different states of 
the system of small particles. An appropriate indicator 
of melting of the fee solid can be defined as 
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FIG. 6: (Color online) Distributions of small particles over 
the rotational invariant we for different scan numbers (indi- 
cated in the figure). The corresponding cumulants We{x) are 
plotted in the upper part of the figure. The position of the 
half- height of W(i{x) can be used as a melting indicator, as 
discussed in the text. 



the value of Wg'^ for the fee lattice. Figure [Sljc) shows the 



values oiM calculated for each scan number. In Ref. [24 1 
we used the threshold value A4 ~ 1.3 at melting. For 
this value, however, both the systems of small and large 
particles appear to melt upon increase in the pressure 
[see Fig. Eljc)]. If we require that only the system of 
small particles exhibits melting (in agreement with the 
other two phase change indicators discussed above) , then 
the "real" melting threshold value should be somewhat 
increased (to A^ ~ 1.5 — 2.0). Qualitatively, the M- 
based measure adequately describes the decay of ordering 
accompanying an increase in the neutral gas pressure. 



THEORETICAL INTERPRETATION 
Plasma parameters and particle charges 

We attribute the observed fluid-solid phase change in 
the system of small particles to the variation in the 
electrical repulsion between the particles. Manipulating 
the gas pressure experimentally, changes various com- 
plex plasma parameters and modifies the strength of the 
repulsion. When the electrical coupling reach (or drops 
below) certain level, freezing (or melting) occurs. To ver- 
ify this scenario we need to estimate relevant complex 
plasma parameters. 



where w^^ is the position of the half- height of Wq{wq) ber 
[so that We{w^^) = 1/2] and u;|'='= = -1.3161 x 10"^ is 



We use the results from SIGL0-2D simulations [25 1 
to estimate plasma parameters in the absence of parti- 
cles. Earlier, it has been shown that SIGLO simulation 
yields reasonable agreement with the results from Lang- 
muir probe measurements for a similar discharge cham- 



56[. In the considered regime (p ~ 10 — 25 Pa, 
rf-amplitude '^ 15 V) the central plasma density is linear 




theory 6ll. |6 
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FIG. 7: (Color online) Estimated values of the particle charge 
number \Q/e\ (a) and effective plasma screening length A (b) 
vs. the scan number. Red circles (blue triangles) correspond 
to the system of large (small) particles. 



on p and can be, with a reasonable accuracy, described 
as TT-o — (1.20 + O.llp) X 10^, where uq is in cm~'^ and 
p in Pa. The electron temperature exhibits almost no 
dependence on pressure, T^, ~ 3.8 eV. Ions and neutrals 
are at room temperature Ti^n '^ 0.03 eV. 

When particles are injected into the discharge they 
inevitably modify plasma parameters. In the follow- 
ing we assume that inside the particle cloud the elec- 
tron temperature remains unaffected, while the electron 
and ion densities are modified to keep quasineutrality, 
rie + \Q/e\np — Ui, where Q is the particle charge. 
Furthermore, we assume that n^ remains close to the 
particle- free value tiq, while n^ somewhat increases in 
response to perturbations from the particle component. 
Physically, this assumption corresponds to the case of 
static equilibrium, when electron and ion densities sat- 
isfy Boltzmann relation in the local plasma potential 
modified by the presence of charged dust particles. It 
is in reasonable agreement with numerical simulation re- 
sults [571, |58j regarding plasma parameters inside the par- 
ticle cloud (simulations show, however, that inside the 
void region in the center of the discharge, the plasma den- 
sity and electron temperature can be considerably higher 
than those in the particle- free discharge [57l|58[). Note 
that in the opposite limiting case, when the ion equi- 
libration length is much longer than the characteristic 
size of the particle cloud (e.g. sufficiently thin cloud), it 
would be more reasonable to assume that the ion den- 
sity is constant, whilst electron density is depleted [59j. 
Such a situation has been apparently realized in a recent 
experiment [6fl] at a very low neutral gas pressure. 

Using these assumption we can calculate the depen- 
dence of the particle charge on pressure (and interpar- 
ticle distance) employing certain model for the electron 
and ion fluxes absorbed on the particle surface and re- 
quiring that these fluxes balance each other. For the 
electron flux we use the orbital motion limited (OML) 



where z ~ \Q\e/aTf^ is the reduced particle charge and 
vTc = yT^JrOe is the electron thermal velocity. For the 
ions, however, the OML approximation would be much 
less appropriate. There are two main reasons for that. 
First, as has been evidenced in a number of recent stud- 
ies [l|, I63l - l7l| , ion- neutral collisions can represent a very 



important factor affecting and regulating the ion flux col- 
lected by the particle. A collision between an ion and a 
neutral atom lowers the ion energy (in this respect reso- 
nant charge exchange collisions are especially important) 
and destroys its angular momentum. If such a collision 
occurs sufficiently close to the negatively charged parti- 
cle, an ion which experienced a collisions has a very low 
probability to overcome the attraction from the particle 
and will eventually reach its surface. Thus, in the con- 
sidered weakly collisional regime, ion- neutral collisions 
result in a more effective ion collection by the particles, 
i.e. they increase the ion flux. The second factor, dis- 
cussed very recently in the context of particle charging in 
plasmas, is associated with the ionization processes. Sim- 
ilarly to ion- neutral collisions, electron impact ionization 
events create slow ions in the vicinity of the particle and 
therefore enhance the ion flux towards its surface [7J]. 

To make a quantitative estimate we employ an expres- 
sion of Lampe et al. [6J| for the collision-enhanced ion 
flux 
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where VTi = yjTi/mi is the ion thermal velocity, r = 
Te/Ti is the electron-to-ion temperature ratio, li is the 
ion mean free path with respect to collisions with neu- 
trals {ti ~ Tn/pain, where a^ ~ 2 x lO"^** cm^ is the 
effective ion- neutral collisions cross section in argon) , and 
i?o determines the radius of a sphere around the particle, 
inside which the potential energy of ion-particle inter- 
action is sufficiently high (higher than the average ki- 
netic energy of an ion after a collision). The first two 
terms in the square brackets on the right-hand side of 
Eq. dH) reproduce the result of the collisionless OML 
model. The third term represents the collision-driven 
flux enhancement. There is no consensus presently re- 
garding the precise determination of Rq in this colli- 
sional term and different arguments have been put for- 
ward 6J,|69|,l7ll|7^. In general, Rq should depend on the 
exact shape of the ion-particle interaction potential, de- 
tails of the ion-neutral collisional proces ses, distribution 
of ion velocities, etc. In a recent paper [7l| it has been 
suggested that in the regime of sufficiently small parti- 
cles and weak ion-particle coupling, the Coulomb radius 
Re — \Q\e/Ti — azT can be used as a relevant measure 
of Rq. This choice has been shown to agree reasonably 
well with the experimental results for very small (1.31 



/xm in diameter) particles [7l|. In the present experi- 
ment, however, the particles sizes are somewhat larger, 
and the applicability condition to set Rq equal to Rq 
(which requires i?c ^ A), is violated, at least for larger 
particles. To get an idea what kind of modification is 
required, let us consider the opposite limit of strong ion- 
particle coupling, Re S> A. In this case, assuming the 
Debye-Hiickel (Yukawa) form of the interaction potential, 
we can easily estimate the length scale of the region of 
strong ion-particle interaction as Rq ~ Aln(i?c/A) [73|. 
A simple expression of the type Rq ~ Aln(l -I- i?c/A) 
would therefore describe adequately the corresponding 
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limits of weak and strong ion-particle coupling and pro- 
vide a smooth transition between them. We adopt this 
heuristic approximation in the analysis below. 

Regarding the effect of ionization, in Ref. [7^ it has 
been demonstrated that the relative magnitude of ion- 
ization and collisional enhancements of the ion flux is 
roughly given by the ratio of the corresponding frequen- 
cies, t'l/fin, where vi is the ionization frequency and 
^'iii ~ VTi/ii is the frequency of ion- neutral collisions. 
This implies that the effects can be added in a simple 
superposition, which yields the following expression for 
the ion flux 



J, ~ V8^a^n,VT, [1 + zr + (1 + lyi/i^i^) (X^/a^e,) W (1 + i?c/A)] 

I 



(10) 



For a fixed gas and ion temperature, the ratio lA/i'in is 
a function of a single parameter - electron temperature 
Tg. The corresponding function has been evaluated for 
neon and argon plasmas with room temperature ions [74 1 . 



7J] we conclude 



Prom the results shown in Fig. 1 of Ref. 
that i^i/fin ~ 1.8 at T^ ~ 3.8 eV. 

The charge can then be estimated from the balance 
condition Ji = J^ using expressions ([8]) and (fTO|) . In 
doing so we take into account the modifications of 
the ion density compared to the particle-free value dis- 
cussed above. Quasineutrality condition implies nijn^ — 
ni/riQ ~ 1 + zP, where P — {aT^/ e'^){np/ tiq) is the scaled 
particle-to-plasma density ratio (the so-called Havnes pa- 
rameter). The effective screening length (screening is 
mostly associated with the ion component for r ^ 1) 
is A = Aq/vT+zP, where Ao = \/Ti/Ane^nQ is the un- 
perturbed ion Debye radius. 

The resulting dependence of the particle charge and 
plasma screening length on the scan number is shown in 
Fig. [T] We observe that the absolute magnitude of the 
particle charge (circles) slightly decreases with increas- 
ing pressure while the plasma screening length (trian- 
gles) exhibits the opposite behavior. The first tendency 
is expected: Increase in the coUisionality lowers the ab- 
solute magnitude of the particle charge in the weakly 
collisional regime [ij, |70|. However, the effect is much 
less pronounced than it would be for an individual par- 
ticle. This is the result of the coupling between Ue, rii, 
and \Q/e\np. Namely, a decrease in Up with increase in 
the pressure results in a decrease in rii and, hence, in 
Ji , so that the collisional enhancement of the ion flux is 
almost compensated by ion depletion in the considered 
regime. The dependence of the screening length on pres- 
sure is unexpected. In the particle-free plasma one would 



expect approximately A ex n, 
of A with increasing pressure. 



-1/2 



the ion density triggered by the expansion of the particle 
cloud with increasing p provides overcompensation and 
the trend reverses. Overall, the relative variations in Q 
and A are rather weak, considerably weaker than those 
in the interparticle distance. 

Thus, the coupling between Ue, n^, and \Q/e\np via 
the quasineutrality condition and the the correspon ding 
effect of charge reduction in dense particle clouds [59|, |60J, 
|75[ plays essential role in the present experiment. Note 
that this effect is to some extent similar to a reduction 
of colloidal charge when increasing the colloidal volume 
fraction. It has been demonstrated that as a result of 
this charge reduction colloids can exhibit an intriguing 
reentrant melting behavior, when a colloidal fluid phase 
appears at a higher volume fraction than a colloidal crys- 
tal 7g,|77|. Similar behavior - reentrant fluid-solid-fluid 
series of phase changes upon isothermal compression - 
can be expected in complex plasmas, as discussed re- 
cently in Ref. [73|- However, the particle clouds observed 
in the present experiment are not dense enough to ex- 
hibit such a behavior. As we have seen, in the regime 
investigated, some increase in A with pressure is accom- 
panied by almost constant values of the particle charge 
and screening length. 



Interaction and coupling strength 

The pairwise potential of electrical interaction between 
highly charged particles is often assumed to be of Debye- 
Hiickel (Yukawa) form, 



C/(r) = (QVOexp(-r/A), 



(11) 
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(X p "'', I.e. a decrease 
Here again, a decrease in 



where r is the interparticle distance. The actual inter- 
actions are known to be considerably more complicated 
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FIG. 8: (Color online) Estimated interparticle interaction en- 
ergy (7(A) at each scan. Blue triangles correspond to a com- 
plex plasma composed of small particles. Red circles are for 
the system of large particles. The dashed curve shows the 
evolution of the neutral gas pressure in the experiment. Note 
that the interaction energy is more than one order of magni- 
tude higher than the gas (room) temperature (r„ ~ 0.03 eV), 
indicating that both systems always remains in the strongly 
coupled state. 



especially at large interparticle separations. In partic- 
ular, continuous plasma absorption on the particle sur- 
face can give rise to unscreened inverse-power-law long- 
range asymptotes of the potential [70, 1791483 1 . In addi- 
tion, plasma openness, associated with constant plasma 
absorption on the particles, can give rise to the so-called 
"ion-shadowing" attraction [79j, i84h87J . It basically rep- 
resents the plasma drag that one particle experiences as a 
consequence of the plasma flux directed to another neigh- 
boring particle and vice versa. Electron and ion pro- 
duction (ionization) and loss (e.g. recombination) can 
result in the emergence of the two dominating exponen- 
tially screened asymptotes (a double- Yukawa repulsive 
potential), one of which is determined by the classical 
mechanism of Debye-Hiickel screening, while the other 
is merely controlled by the balance between the plasma 
production and loss 88|, |89| . Here we neglect all these 
corrections and adopt the simple form (|11|) . The justify- 
ing arguments are as follows: (i) The mean interparticle 
separation is not large enough in the present experiment 
(A ^ 3 A) for the IPL asymptotes to dominate [82J; (ii) 
The neutral gas pressures used in this experiment are well 
above those required to make "ion shadowing" attraction 
operational '9u\; (iii) The relative ionization efficiency is 
not high enough to expect significant deviations from the 
conventional screening regime (lll[) at distances charac- 
terizing the mean interparticle separation [891]. 

Using the measured values of A and estimated values 
of Q and A we can evaluate variations in the interaction 
energy between neighboring particles U{A). The results 
are shown in Fig. [5] The energy, as a function of p, ex- 



hibits maximum at the lowest pressure p ~ 11 Pa in both 
cases studied. Thus, increasing (decreasing) the neutral 
gas pressure implies reduction (grows) in the strength of 
electrical repulsion (coupling strength) between the par- 
ticles. As discussed above, the main factor responsible for 
this behavior is the compression (expansion) of the par- 
ticle system upon decrease (increase) in the pressure (Q 
and A are less sensitive to pressure manipulations). Con- 
sequently, complex plasma appears more ordered when 
the pressure is low and less ordered when the pressure is 
high (in the pressure range investigated). Freezing and 
melting can also be expected, provided the system is lo- 
cated not too far from the phase boundary. This is in full 
qualitative agreement with the results from the structural 
analysis. 



Equilibrium Phase Diagram 

To get further insight into phase behavior of complex 
plasmas under investigation, let us estimate the location 
of phase trajectories of these system on a relevant equi- 
librium phase diagram. 

The system of particles interacting via the Yukawa 
potential (jlip in thermodynamical equilibrium can be 
characterized by two dimensionless parameters. In the 
field of complex (dusty) plasmas the screening param- 
eter K — A/A, characterizing the efficiency of screen- 
ing, and the coupling parameter T = (Q^/TpA) exp(— k), 
measuring the ratio of the interaction energy at the mean 
interparticle separation to the particle kinetic tempera- 
ture, are commonly used. Equilibrium phase diagrams 



of Yukawa systems have been extensively studied [91- 
|95| . In the strongly coupled regime they can exist in the 
fluid or solid phases. In the solid phase particle form ei- 
ther bcc (weak screening regime) or fee (strong screening 
regime) lattices. The triple point is located at k ~ 6.9 
and r ~ 3.5 [95| (according to a more recent estimate 
from Ref. [96| its location is at k ~ 7.7 and F ~ 3.1). The 
boundary between the fluid and solid phases on the plane 
(k, F) can be approximated by the expression [97|, |98j 
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1 



1^2- 



(12) 



where the subscript "M" refers to melting (the density 
gap is rather small so that it makes little sense to dis- 
tinguish between freezing and melting here). Other ap- 
proximate analytical expressions have been also proposed 
in the literature to locate the fluid-solid coexi stence for 
a wide class of interaction potentials 99|-[l01[. We use 
Eq. ([12]) here due to its particular simplicity and rea- 



sonable accuracy: Deviations between the results from 
Eq. (fT2|) and numerical simulation data from Ref. 95| do 
not exceed several percent, as long as k remains not too 
large (k < 8) ^. 
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FIG. 9: (Color online) Estimated trajectories of complex plas- 
mas, studied in the present experiment, on the equilibrium 
phase diagram of Yukawa systems. The solid curve corre- 
sponds to the fluid-solid phase change [Eq. (I12|) ]. the dashed 
curve shows the boundary between the bcc and fee solids 
(smooth fit to the simulation data points from Ref. [93] )• Blue 
triangles denote phase state points for the system of small 
particles, red circles correspond to large particles. Arrows 
mark the direction of the phase evolution when the pressure 
increases. 



Figure [S] shows the equilibrium phase diagram of 
Yukawa systems in the (k, F) plane along with the es- 
timated phase-state points from the present experiment. 
In this figure circles and triangles correspond to the phase 
states visited by the system of small and large particles, 
respectively. The arrows indicate the direction of phase 
evolution when the pressure increases. With increasing 
pressure both systems move towards the melting curve. 
Complex plasma composed of small particles crosses the 
phase boundary for highest pressures investigated. At 
the same time, complex plasma formed by large parti- 
cles always remains in the solid state. These result are 
in reasonable quantitative agreement with those from 
structural order measurements and estimated values of 
freezing/melting indicators. On the other hand, there 
is a notable disagreement with respect to the structure 
of the solid phase. Equilibrium Yukawa crystals, in the 
range of n investigated, should form the bcc lattice (see 
Fig. ini). In contrast, the observed crystalline structures 
are dominated by the hep and fee lattices (see Fig. 2]) . 
Such a drastic discrepancy can be associated with the 
fact the crystallization evolves over several essentially 
non-equilibrium stages, yet the physics of these non- 
equilibrium processes is still poorly understood. In par- 
ticular, this concerns the crystal nucleation - both homo- 
geneous and heterogeneous - out of the supercooled melt 
(when the thermodynamically stable solid has a s truct ure 
which is vastly different from that of the liquid) |102| , as 
well as the crystal nucleation out of a thermodynamically 



unstable solid phase (such as hep and fee lattices, which 
have a structure incompatible with the ground-state bcc 
crystal). In this case, one could expect largely prevented 
nucleation of the equilibrium phase, and our experiment 
might be a manifestation of this effect. 



DISCUSSION 

In the first publication related to the discussed exper- 
iment [2^ the theoretical interpretation was somewhat 
different. In particular, we did not take into account the 
effect of ionization enhanced ion collection by the par- 
ticle. In this way, the particle charge was apparently 
overestimated. As a result, we had to assume that the 
particle kinetic temperature Tp was several times higher 
than the neutral gas temperature in order to explain 
melting. Estimates presented in the present paper yield 
lower charges and, therefore, make such an assumption 
unnecessary (see Fig.[9|). Unfortunately, the particle dy- 
namics is not resolved in this experiment, so we cannot 
completely exclude the possibility that Tp is somewhat 
higher than r„. 

It is instructive to discuss some specific properties of 
the phase transitions observed. As pointed out above, 
an increase in A (i.e. decrease in the particle density) 
is the main factor responsible for the melting when the 
neutral gas pressure increases. Similarly, reducing pres- 
sure compresses the particle system and this stimulates 
freezing. This is a generic mechanism of (isothermal) 
fluid-solid phase transition, that can be realized in ex- 
tremely wide range of various substances and materials. 
In complex plasmas, it can be in principle observed in 
both 3D and 2D systems 18|, |22|. However, 2D complex 
plasmas are known to be strongly affected by plasma- 
specific mechanisms of melting. One of the clear man- 
ifestations is the conventional procedure of melting flat 
plasma crystals by reduci ng the n eutral gas pressure in 
ground-based experiments [13|,|l03[. The difference is not 
a consequence of the essentially 2D character of crystals 
investigated on Earth, but is rather due to the presence 
of strong electric fields (and, therefore, strong ion fiows) 
required to balance the force of gravity. There are ef- 
fective mechanisms of converting the energy associated 
with ion flows into the kinetic energy of the particles. 
Know n scenario s include ion-particle two-stream insta- 
bility [104L |l05|, non-reciprocity of the interaction due 
to asymmetric character of the scre ening clou d around 
the particles ( "plasmawakes" ) 2l|, liod llOTf . particle 
charge variations 108l , ll09| . All these scenarios lead to an 
abrupt increase of the particle kinetic energy at pressures 
below certain threshold value, causing crystal melting. 
The process of melting observed in the present experi- 
ment is apparently free of these plasma-related effects, 
but has much more in common with generic processes in 
conventional atomic, molecular, and soft matter systems. 
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CONCLUSIONS 

We have presented a comprehensive overview of the 
experimental studies on fiuid-sohd phase transitions in 
large 3D complex plasma clouds performed in micrograv- 
ity conditions onboard the ISS. The neutral gas pressure 
turns out to be a convenient control parameter to drive 
crystallization and melting. In the parameter regime in- 
vestigated, the phase transition is mostly associated with 
the compression (expansion) of the complex plasma sys- 
tem upon decrease (increase) in pressure, associated with 
the variations in plasma confinement. This is very differ- 
ent from the conventional procedure of melting flat (2D) 
complex plasma crystals by reducing the gas pressure in 
ground-based experiments. Detailed analysis of complex 
plasmas structural properties and evaluation of three dif- 
ferent freezing/melting indicators reveal an overall good 
qualitative agreement. Theoretical estimates of the com- 
plex plasma parameters allow us to approximately de- 
termine the locations of the systems investigated on the 
relevant (Yukawa systems) equilibrium phase diagram. 
Here again, reasonable agreement is documented, except 
the structure of the solid phase, which is expected to 
be bcc lattice form numerical simulations, but is dom- 
inated by the hep and fee domains in experiment. We 
attribute this to the non-equilibrium character of phase 
evolutions studied in the present experiment, but this 
issue certainly deserves more attention. Further investi- 
gations are planned, which we believe will provide impor- 
tant insight regarding the principle mechanisms dominat- 
ing ubiquitous and still poorly understood phenomena of 
crystallization and melting in complex plasmas and re- 
lated systems of strongly coupled particles. 
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